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Bacher and Goudsmit Theory of Complex Spectra 

R. E. Trees 

The Bacher and Goudsmit theory expresses the energy of a state of an atom in terms 
of a linear combination of energies of states of other ions of the same atom. This theory 
is tested in the spectra with d- and s-electrons and general conclusions drawn about its 
accuracy. 

In deriving their theory, Bacher and Goudsmit formally introduced many-particle 
interactions. They showed (to a limited extent) that this was equivalent to a second 
procedure in which the zero-order wave functions were assumed to be independent of the 
ionization and successive orders of perturbations were included. 

Recently it has been shown that a linear theory, which is a combination of the Slater 
theory and an empirically derived correction, predicts the separations of terms in any one 
configuration with about the same accuracy as the Bacher and Goudsmit theory; in this 
theory, the wave functions are allowed to depend on the degree of ionization. It is shown 
that the Bacher and Goudsmit theory can be interpreted in a third way that is consistent 
with this linear theory. 

In the light of this latter interpretation, Racah has si invested an empirical modification 
of the Bacher and Goudsmit theory that often leads to better results; this suggestion has 
been incorporated in the present paper. 



1. Introduction 

Bacher and Goudsmit have shown that it is pos- 
sible to express the energy levels of an atom cor- 
rectly to second order by taking linear combinations 
of the energy levels of the ions of that atom. 1 They 
gave explicit formulas for these linear combinations 
for atoms containing s- and p-electrons, and tested 
the theory in the spectra of oxygen, nitrogen, and 
carbon. The theory usualty predicted the absolute 
term values with an error less than 2500 K. 2 Sepa- 
rations of terms in a given configuration were often 
given very much better, sometimes within a few 
hundred kaysers. 

In section 2 the theory of Bacher and Goudsmit is 
derived briefly from a new viewpoint. Bacher and 
Goudsmit derived their formulas by formally intro- 
ducing many-body interactions; no explicit form was 
given for these interactions. They were able to 
show that in the first and second approximations this 
was equivalent to using zero-order wave functions 
that were independent of the degree of ionization, 
and, with these, making first- and second-order per- 
turbation calculations of the energy, respectively. 
In the present paper, a third derivation is given, in 
which the zero-order wave functions of the Slater 
theory are allowed to vary regularly with the degree 
of ionization; an explicit two-body interaction is also 
introduced. This interpretation is justified by the 
observation that there is a smooth variation of the 
one-particle and two-particle parameters in a linear 
theory. 3 The more specific and easily understood 
meaning of this interpretation in terms of a linear 

i R. F. Bacher and S. Goudsmit, Phys. Rev. 46, 948 (1934); referred to in the 
present paper as "BG". 

2 "Kayser" is the name adopted by the Joint Commission for Spectroscopy for 
the unit of wave number, cm-', and is abbreviated "K". 

3 R. K. Trees and M. M. Harvey, J. Research NBS 49, 397 (1952) RP2378. In 
section 1 of R P2378 the literature is reviewed, and the methods of analyzing spec- 
tra with the Slater theory, for including effects of configuration interaction, and 
of making allowance for polarization energy with the L{L-\-l) correction are also 
outlined and applied. In referring to a "linear theory", the L{L-\-\) correction 
is combined with the Slater theory, and the radial integrals are regarded as 
parameters. 
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theory is given in section 5; it is equivalent to assum- 
ing that the parameters in different ions of the same 
atom are given by series expansions in powers of the 
degree of ionization. It is shown that the BG nth 
approximation accounts for terms up to and including 
the nth power for the one-particle parameters and 
terms to the (n— l)th power for the two-particle 
parameters. 

In section 3 a few formulas applicable to terms in 
d n , d n s, and d n s 2 configurations are obtained by using 
tables of coefficients of fractional parentage ami 
recoupling transformations given by Racah. 4 In 
section 4 an application of these formulas is made in 
the spectra of titanium, vanadium, chromium, and 
zirconium. It is found that separations of terms in 
the same configuration are given very closely (the 
error is less than 400 K). Separations of terms in 
different configurations may be given very poorly, 
however, as an error of 6000 K is found hi titanium 
(Z=22). The latter type of error is found to 
decrease very rapidly in spectra of more highly 
ionized atoms: in chromium spectra (Z=24) the 
error is only 600 K, one-tenth the error found for the 
related titanium spectra. The error is also found to 
decrease for spectra with higher nuclear charge as it 
is only 2000 K in zirconium (Z=40). 

Racah has pointed out that a power-series expan- 
sion, of the type described previously, converges 
more rapidly if spectra are chosen having the same 
degree of ionization. A comparison with experiment, 
using this empirical modification of the original 
theory of Bacher and Goudsmit, is also given in 
section 4. 

As indicated previously, the Bacher and Goudsmit 
theory can be interpreted in terms of a linear theory 
in which the radial parameters are given by series 
expansions in the degree of ionization. Checking 
these expansions is the easiest way to check all the 
data against the Bacher and Goudsmit theory. The 



« G. Racah, Phys. Rev. 62, 438 (1942); 63, 367 (1943). These papers aie referred 
to herein as II and III, respectively. 
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conclusions reached from the direct application of the 
theory already cited above, have been verified by 
using all available data with this method. The de- 
tailed calculations needed for this verification are 
omitted for brevity. 

Jn section 6 the effect of weak configuration inter- 
action in which one electron is excited is considered, 
and in section 7 the effect of irrationality of the 
energy eigenvalues is analyzed. The purpose is to 
illustrate in more detail the direct application of the 
theory of BG and to indicate where a direct applica- 
tion might lead to results that differ from those 
obtainable with a linear theory. Such differences are 
probably too small to be detected within the accuracy 
obtainable at present with any theory, but more com- 
plete comparisons of theory with the presently avail- 
able experimental data are needed to verify this 
conclusion. 



2. Bacher and Goudsmit Theory 

2.1. General Theory 

In this section is given a new derivation of the 
theory in which it is assumed that the interaction 
parameters vary smoothly. The Hamiltonian of 
the 7i-electron system is taken as 



where 
and 



H=H 1 + H 2 , 



(1) 



*=s(£*-f) 



i>j = l 'ij 



Hi and H 2 are examples of what are called " type- 
one/ 7 and "type-two" terms, respectively, {H x is a 
sum of functions, each depending on the coordinates 
of a single electron, and H 2 is a sum of functions for 
two electrons) . A particular selection P of r of the 
^-electrons will be indicated by P r . The sum of the 
Hamilton ians H(P r ) of the r-electron systems for all 
possible choices P is a symmetric function of the 
^-electrons, and the ^-electron Hamiltonian H is 
expressible identically as a linear combination of such 
symmetric sums for two or more values of r, 






(2) 



The conditions that this be an identity are the follow- 
ing: 

?5CXD=G) ^ =<)©-G> < 3 > 

In general, ( , ) is the binomial coefficient when 
a>&>0, 



fa\_ a\ 

\b)~{a-b)\ 



b\' 



and zero for 6<0 or a<Cb. The values of the co- 
efficients C T will be indeterminate if there are more 
than two nonzero (7 r -values. The method of choos- 
ing best values of the C r in such cases will be con- 
sidered below. 

The eigenfunction ^(7) of the ^-electron system is 
assumed to have in its specification a definite con- 
figuration l h l 2 , . . . } l n , and P T is specified by r of these 
Z-values; the electrons are regarded as distinguishable, 
even though some may be equivalent. The electrons 
in the selection P r have allowed eigenf unctions \p(a) 
of the r-electron Hamiltonians H(P r ); the residual 
(n—r) electrons have allowed eigenf unctions \p(fi) of 
the (n— r)-electron Hamiltonian. The /-values of 
the electrons are included in the specifications a, /?, 
and 7, along with S- and L- values for the combined 
system and any other quantum numbers needed to 
specify the coupling. 5 The function formed by vec- 
tor coupling the (n— r) -electron eigenfunction to 
the r-electron eigenfunction is indicated by \j/(a,l3,y). 
A linear combination of thes? functions can be formed 
that will be identical with ^(7) ; this was first pointed 
out by Bacher and Goudsmit, 



Hy)-- 



= S(«/MhO*(«,/3,7). 

a/3 



(4) 



(a,]3,7}7) is called the coefficient of fractional parent- 
age. The calculation of the coefficients of fractional 
parentage is a mathematical problem that has been 
considered in detail by Racah. He has tabulated 
these coefficients for configurations of equivalent p- 
and d-electrons in III; their use will be illustrated in 
section 3. 

A partial energy E(a) of the choice P r is introduced 
that is associated with the r-electrons while they are 
part of the ^-electron system and combined to form 
the state a. This energy cannot be observed, and 
must be approximated from the observed energy of 
the r-electron ion of the atom in the state a, which is 
defined as W(a). The average energy of the choice 
P r when present in the original 71-electron configura- 
tion 7 is determined with the help of (4), 



a/3 



(5) 



E(y)=TiCrT;E(r,Pr)- 



(7) 



V 



5 There may be more than one term allowed, with a given (SL) -value in the 
configuration considered (the configuration may be that contained in a, /S or y). 
In such cases the function obtained by a legitimate coupling procedure will not 
be an eigenfuncticn of the Hamiltonian, even though it has an allowed symmetry, 
if it mixes with other allowed eigenf unctions in that configuration having the 
same (SL) -value. In such cases a purely formal difficulty arises because the 
characterization by parentage is indeterminate, and the theory given by Bacher 
and Goudsmit must be amplified by more or less arbitrary additions. In sec- 
tions 2, 3, and 4 of the present paper these cases are omitted. These omitted 
cases are considered in sections 5 and 7. 



The corresponding approximate energy can be cal- 
culated from the observed energies in the r-electron 
ion, 

W{r,P r )=12(etPnh) 2 W{a). (6) 

or/3 

The energy of the ^-electron system is then given 
identically by an expression similar to (2), 
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The corresponding approximate energy is calculated 
from' the approximate average energies (6), 



r v 



(8) 



The energies W and W used above are the same as 
the W used in BG for both expressions. Quantities 
such as W(a), W(a,P), etc., on the right side of eq (7) 
to (14) in BG are sometimes energies of observed 
spectral terms, but in general they are linear com- 
binations of the observed energies; in the present 
notation they would be replaced by W(y;Pi), W(y;P 2 ), 
etc. Directly observed energies are indicated by W, 
as in the right-hand side of (25) and (26) of BG. 

The difference between the exact value (7) and the 
approximate value (8) is estimated by separating 
the E's and W's into their type-one and type-two 
parts, and by assuming that the following expansions 
can be made: 



W 1 (a)=E 1 (a){l+a 1 (y)(n-r) + a 2 (y)(n-ry+. . .] 



(9) 



W 2 (a)=E 2 (a){l + b 1 (y)(n-?')+b 2 (y)(n-r) 2 + . . .}. 

The coefficient a t (y) in (9) will be eliminated from 
the difference between (7) and (8) if 

g^- r K")(0 =o ' (io) 

and the coefficient b t (y) will be eliminated if 

§«*-'>• C)(a)-o. ("> 

A choice of C r - values defines an approximation of 
the same order as the largest value of t for which 
(10) is satisfied (it is assumed that (10) and (11) are 
satisfied for all smaller values of t and that relations 
(3) are satisfied) . For instance, in a second approxi- 
mation, the type-one energies are expanded in a 
series containing powers up to and including the 
squares of (n—r), whereas the type-two energies are 
expanded only to first power in {n—r). Formulas 
(3), (10), and (11) are not independent. They may 
be combined into a definition that is easily shown to 
be equivalent to that just given, by requiring the 
coefficients C r to satisfy the following (m+1) equa- 
tions in order to define an mth approximation 



§**-<)©-"© 



0=0,1,..., m). 

(12) 



In the derivation of BG, many-body interactions 
are formally introduced. According to the present 
notation, this is equivalent to considering not only 
type-one and type-two terms as being present in (1) 
but also type-three, type-four, etc., terms. This 
leads immediately, by generalizing (3), to the defini- 
tion of the m th approximation given by Bacher and 



Goudsmit. According to their definition the re- 
values satisfy the following (m+1) equations: 

§<)GD = © t = 1 '2,..,(m+l). (13) 

To show the definitions (12) and (13) are equivalent, 
the following expansion is utilized for powers £>0 of 
an integer p: 6 



where 



a.®- 



p'=i>«(o(f)> 



(14) 



Inserting (14) in (12) leads to the result 

»^<: 1 )g<:r 1 )-« *-'• 

Equation (13) is equivalent to the formula 



.,m). 
(12a) 



s<=*> 



= 1 (£=1,2,..., (m+1)). (13a) 



It follows therefore by repeated application of the 
binominal coefficient relation, 

/n — k— l\_/w — k\ hi — k— 1\ 
\r-l J~\r-l)\r-2 ), 

and the use of (13a) that the left side of (12a) 
vanishes, as it should, whenever (13) is satisfied. 
When 2=0 in (12), the equation is the same as (13) 
when k=l, so the equivalence of the two sets of 
defining equations is demonstrated. A special m th 
approximation is defined by the following set of C T - 
values: 

C,=(-)-' +1 (^~j) [fr-1,2,... ,<«+!)) 



a=o 



(r=(m+2) > ...,(»-l))J 

(15) 



If there are only two nonzero values of C r , and 
these are ft and C 2 , then they are defined by (15), 
and relation (8) has the explicit form 



Wi(7) = S W(r,P2)-(n- 2E^( 7 ;Pi). (16) 

V V 

Relations (8) and (12) of BG are special cases of 
(16), for which n has the values 4 and 5, respectively. 
As an alternative, it can be assumed that d and 
C n -\ are the nonzero coefficients, and then the fol- 

6 A treatise on the calculus of finite differences, G. Boole, edited by F. J, 
Moulton (Macmillan & Co., London, 1880). 
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lowing formula is obtained by solving (12) or (13) 
(the specification y will be omitted for brevity): 



W 1 (n) 



1 ^VTTV 



n z p 



^SWO- (17) 



2.2. The Second Approximation 

If Ci, C 2 , and C 3 are ta^en as the nonzero coeffi- 
cients of the second approximation, the following 
formula is obtained from (15): 



W 2 (n)= 2 W (P 3 ) - (k - 3)1] W(P 2 ) + 



1 



(n— 3)(n- 



•2)SW'(Pi). 

P 



(U 



Relations (9) and (13) of BG are special cases of this 
formula for which n has the values 4 and 5, respec- 
tively. In appendix I of BG, perturbation theory 
was used to show that (18) is equivalent to a calcu- 
lation that is correct up to and including the second 
order. In general, reference to "the theory of 
Bacher and Goudsmit" in the present work applies 
chiefly to this second approximation. 

The complete set of one-electron functions used 
in appendix I of BG was that appropriate to the 
ionized state of the atom with only a single electron 
outside of a closed shell. The same zero-order 
functions were used for the atom in all other states 
of ionization, and the second-order effect was evalu- 
ated with them. This is a poor zero-order approxi- 
mation for some spectra, 7 so that "second-order 
effect," in the sense used by Bacher and Goudsmit, 
is not synonymous with polarization energy. The 
latter energy is usually defined as that part of the 
observed energy that cannot be explained by a cal- 
culation in which the one-electron functions are 
determined by the Hartree-Fock method for each 
state of ionization. No investigation has been made 
to determine the variation in accuracy of a second- 
order calculation with change of the zero-order 
eigenf unctions, but Bacher and Goudsmit often 
obtained poor values in their calculations because 
the higher-order effects became more important as 
the zero-order functions became worse approxi- 
mations. 

The influence of polarization on the separation of 
terms in a single configuration containing d- and 
s-electrons may be represented by a correction pro- 
portional to L(X+1). 8 Racah has shown that this 
correction may be derived by assuming a special 
type of linearity for the polarization energy, and he 
bias outlined a proof of this assumption of linearity. 9 
The proof is based on the assumption that in the 
configuration d n the L(L+1) correction arises from 

7 Similarly defined parameters change by 10 to 30 percent in successive ions 
when they are evaluated by least squares; they would have the same value if the 
zero-order functions were the same in all ions. A similar problem is considered 
by A. Rahman, Physica 19, 377 (1953). See also E. P. Wigner, Phys. Rev. 94, 77 
(1954). 

s R. E. Trees, Phys. Rev. 84, 1089 (1951). 

9 G Racah, Phys. Rev. 85, 381, (1952), and private communications. 



interactions with configurations d n ~ 2 l 2 and d n ~ 2 W that 
have two electrons of the original configuration 
excited. The linear behavior then follows from (54) 
of BG. Alternatively, it can be demonstrated by 
comparing (33a) of III with (33c) of III, and assum- 
ing that the perturbing configurations are suffi- 
ciently far away so that the energy denominators in 
the second-order perturbation can be assumed con- 
stant. The proof shows that the i(Z+l) correction 
is accounted for initially in the first approximation, 
so that it is on the same footing as the Slater energy. 
An interaction in which only one electron is excited 
is included fully for the first time in the second ap- 
proximation, as shown by (53) and (59) of BG. 
Any deviation from linear behavior that is produced 
by such an interaction must be assumed negligible 
if the Z(Z+1) correction is used to account for 
polarization. 

In section 4 it is shown that the Bacher and Goud- 
smit theory leads to close agreement with the ob- 
served separation for two terms in the d* s con- 
figuration in both titanium and zirconium, the errors 
being 357 K and 293 K, respectively. Bacher and 
Goudsmit noted that the separation of terms in the 
same configuration are well predicted in spectra 
containing 2s- and 2p-electrons. 10 It is concluded 
that the error of the Bacher and Goudsmit theory in 
predicting separations of terms in the same con- 
figuration will generally be less than 400 K and that 
the accuracy is therefore about the same as that ob- 
tained with a linear theory. 

A large part of the polarization energy is the same 
for all terms of a given configuration. The energy 
common to all terms can be investigated with self- 
consistent field and other variational methods, but 
it is not easy to get a reliable picture of the behavior 
from the few calculations that have been made. 
One approximation, used fairly often with minor 
variations, ascribes an energy of 5000 to 10000 K 
either to each electron or to each pair of equivalent 
electrons; interactions between pairs of nonequiva- 
lent electrons are neglected in the latter approxima- 
tion. The energy is independent of Z but depends 
a little on the types of electrons that are paired. 11 
This part of the polarization energy is very large, 
and the departures from the described behavior 
could be appreciable. 

The Bacher and Goudsmit theory can be tested in 
this respect by calculating the separations of terms 
in different configurations. In section 4 it is shown 



i° Tables XVII and XVIII of BG, calculation II. It is noted in BG that 
agreement better than 200 to 300 K must be regarded as fortuitous so that the 
ultimate significant accuracy is regarded as about the same as obtained with the 
L(L-j-l) correction. Many separations of terms in the same configuration are 
predicted with errors greater than 1000 K in BG, and this must be attributed to 
the neglect of configuration interaction, or to the fact that the theory is less ac- 
curate for spectra with 2s- and 2p-electrons. 

ii F. Seitz, Modern theory of solids, p. 247ff (McGraw-Hill Book Co., New 
York, N. Y., 1940). Seitz evaluates energy per electron rather than per pair of 
equivalent electrons, but his remarks on page 248 justify the latter interpretation. 
Self-consistent field calculations of O, 0+, and ++ , including exchange and con- 
figuration interaction (D. R. Hartree, W. Hartree, B. Swirles, Phil. Trans. Roy. 
Soc. [A] 238, 229 (1939)), support an interpretation based on pairs, a value of 7000 
K being associated with each pair of equivalent 2p-electrons. Although a strong 
dependence on multiplicity is expected (page 250), this is shown only to a minor 
extent in the oxygen calculations; if present, it would be accounted for in the 
previous consideration of the separations of terms in the same configuration. An 
additional point that should be noted is that comparison of the self-consistent 
field calculations in oxygen (noted above), with and without configuration inter- 
action, indicates that configuration interaction is important. 
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that there is an error of 6000 K in predicting the 
energy difference between a term in the 3d} con- 
figuration and one in the 3d 3 4s configuration of tita- 
nium, and a corresponding error of 2000 K in similar 
zirconium spectrum with 4d- and 5s-electrons, and 
also that this error decreases rapidly as ionization 
increases in the spectra with 3d- and 4s-electrons 
because the spectra of titanium, vanadium, and 
chromium show errors of 6000, 1800, and 600 K, 
respectively. Bacher and Goudsmit found an error 
of about 2000 K in the spectra with 2s- and 2p- 
electrons; the decrease of the error with increase of 
Z was slight in these spectra. 12 

The spectra considered here contain four electrons 
outside of closed shells and are the easiest spectra to 
analyze with relation (18). A more complete ex- 
ploration of the theory would consider also spectra 
with five or more electrons. 

2.3. Higher Approximations 

A third approximation may be obtained from (15), 
in which d, C 2 , (7 3 , and C 4 are nonzero; 

W,(n) =2 W(P t )- (« - 4)2 W(P Z ) + 



<"-4)(n-8) aF(Pi) _ 



(ft-4)(n-3)Q-2) y^WiP,). 
o . 2 n 



(19) 



Relation (14) of BG is a special case of this formula 
for which n=d. Bacher and Goudsmit considered 
this approximation equivalent to a calculation by 
perturbation theory correct to thrd order, but they 
were not able to carry out the complex calculation 
that would be necessary to demonstrate this rigor- 
ously. Separations of terms in spectra of oxygen 
and nitrogen were calculated a little better in the 
W 3 than in the W 2 approximation. 13 In section 5 the 
relative accuracy of W z and W 2 is discussed from 
another viewpoint (formulas (56) and (60)), and it is 
concluded that higher approximations can only ac- 
count for a small part of the error already noted in 
the W 2 approximation. 

Any interpretation used to justify the method of 
Bacher and Goudsmit leads to the conclusion that 
it is more accurate the more similar the binding of 
the electrons (and the structure of the configurations) 
for the different energy levels utilized in the formulas. 

12 Table XIV of BG. In discussing this table, Bacher and Goudsmit noted 
that the errors decreased slightly as 7, increased, but from a later discussion it 
is apparent that they considered this effect of Z as being of minor importance. 

13 Separations of terms in the same configuration were sometimes given with 
i -mailer error by W3 than by Wi, as may be seen by comparing calculation III 
with calculation II in table XVIII of BG; however, in an almost equal number 
of cases the agreement was worse. Absolute term values calculated with Wi for 
the 2.S- 2 2/> 3 configuration of oxygen were generally in error by less than 10C0 K 
(table XLV of BG), whereas the Wz approximation indicated term values in the 
2s 2 2p 4 configuration with a larger error, about 1400 K (table XV of BG) . Bacher 
and Goudsmit note that the 2s 2 2p* configuration is calculated better in the Wi 
approximation, regarding the ,s-clectrons as part of the closed shell, than in the 
Wi approximation. The relative errors of absolute term values will naturally 
decrease for the higher approximation as the absolute value of the energy to the 
closed shell is much greater. 



For example, in calculating term values for l n all 
configurations specifying a choice of /-elections would 
belong to l r . On the other hand, in calculating term 
values for l n ~ l s the selection of /-elections is specified 
by either the l T or the l T ~ l s configuration, so that 
there are twice as many different configurations to 
be considered. According to this argument, a given 
order of approximation would probably predict terms 



in the l n configuration best, terms in 



next best, 



and terms in l n ~~ 2 s 2 least accurately. This argument 
has some confirmation which was pointed out by 
Bacher and Goudsmit. 14 Bacher and Goudsmit also 
have noted that the method is inadequate for con- 
figurations such as 2s 2 2pHs, in which the total energy 
of the excited electron is less than the electrostatic 
energy of the inner ones. In such cases, even the 
separations of terms in the same configurations would 
not be given accurately. Similiarly, the method gives 
poor predictions for negative ions. 15 



3. Formulas for Spectra With el- 
and s-Electrons 

3.1. The 1" Configuration 

When the configuration specified in y is l n , equation 
(()) has the form 

W(/ n ySL;n= , S (l r (a 1 S 1 L 1 )l n ^(a 2 S 2 L 2 ) f 

SL}l n 7 SL) 2 Wil'cuSiLJ. (20) 

If this relation is used in (17) for W(P n -\), the sum 
over P introduces a factor ( j^f ] = /?, and (17) 

becomes, 

W,{l n ySL)=-^- Tj (^-'(a.S.L^ISLjl-ySL) 2 

71 L aj S\Li 



•W^ n " 1 aiS'ii.)- 



n 



n—2 



W(Q. (21) 



The first term on the right defines the separation of 
terms in the l n configuration; it has the same form 
as a relation given by Racah for calculating matrix 
elements of a general scalar operator G defined be- 
tween pairs of interacting electrons. 16 Similiarly, 
by utilizing (20) with r=2, relation (16) becomes 

W^ySL)^^ 2 {IXS.LW-XatStU), 

£ oj ■/>! c*2 ^2 ^2 

SL}l"ySL) 2 W{V l S l L l )-n{n—2) W(l). (22) 



14 Table XIV of B G. Absolute term values for the 2p 4 configuration of oxygen 
are predicted within a few hundred kaysers, which may be the best significant 
agreement obtainable. Errors for 2p 3 2s are about 1000 K, and they are about 
2000 K for 2«2 2p*. 

15 In footnote 7 of BG it is noted that the calculated electron aflinity of oxygen 
is 8000 K and that the observed value is 18000 K. The discontinuity in binding 
energy of negative iens has been discussed by Ta-You Wu in Phys. Rev. 89, 629 
(1952). 

i« G. Racah, Phys. Rev. 76, 1349 (1949). Formula (1). 
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This is the same as an alternative formula given 
by Racah for calculating matrix elements of a scalar 
operator. 17 

It is noted in III that the formula corresponding 
to (22) needs very long calculations and is of practical 
use only in a few particular cases. Although, the 
direct use of (20) is impractical, except when r=n—l 
(or in the trivial case where r is unity) , a step-by-step 
procedure, starting with r=n— 1 and decreasing r 
a unit at a time, is generally practical. With this 
procedure, only tables of the coefficients (l n ~ l 
{a x SiL x )l ,SL}l n ySL) are needed. These coefficients 
are defined by (9) III and (11) III, and the coefficients 
have been tabulated in III for p n and d n configu- 
rations. 

As an example, the second approximation (18) 
will be evaluated in detail for a d A 5 D term by using 
(20) in a step-by-step procedure. By using table 
III of III, relation (20) becomes 

W(d A 5 D;d 3 ) =3/10W(d 3 4 P)+ 7/10W(d 3 4 F) . (23) 

The W's on the right are formally evaluated with 
(20), by using table II of III: 

W(d z 4 P;d 2 ) = 8/15lT(d 2 3 P)+ 7/15W(d 2 3 F), (24) 

W(d 3 ±F;d 2 ) = l/5W(d 2 3 P) + 4/5W(d 2 3 F) . (25) 

These results are then substituted back in (23) : 

W(d 4 5 D;d 2 )=3/10W(d 2 3 P)+7/10TTW 2 3 F). (26) 

Alternately, (26) could have been obtained directly 
from (20) (with r=2) by first calculating the coeffi- 
cients of fractional parentage (d 2 (SiLi)d 2 (S 2 L<>), 
5 Djd' 5 D) with (32) of III. The step-by-step 
procedure used above is simpler because the quantum 
numbers a 2 S 2 L 2 do not enter the calculation explic- 
itly. The values (23) and (26) are now inserted 
into (18) and the sums over P replaced by factors 

( ) to obtain the desired formula, 

W 2 (d* 5 D) = 6i5W(d* 4 P)+ 14/5 W(d d 4 F) 

-9/5 W 2 3 P)-21/5W T (^ 2 3 F)+4W(d). (27) 

Formula (27) is a good one to use in testing the 
Bacher and Goudsmit theory. AH term values used 
in the formula should be unperturbed, and they are 
likely to be known experimentally. The following 
two formulas, obtained in a similar fashion, are 
not so generally useful: the terms that are most 
likely to be perturbed by configuration interaction 

17 Formula (33a) of III. To get the expressions to agree exactly, the following 
generalization cf formula (29) of III is used: 



(WaiSiLi)l"-r(a2S 2 L2)SL]l»ySL) 

= (—l)S+L+L 2 -, 



Si-L l +iS l 8 2 (ln-r( a2 S 2 L2)WaiSiLi),SL]l»ySL) 



are starred : 

W 2 (d* 1 I) = 6/5W r (cP 2 G)*+ 14/5 W (d 3 2 H) 

-3/7 W 2 l T>)*-ZW(d 2 3 F)-l8/7W(d 2 l G) + 4W(d) 

(28) 

W 2 (d" 3 H)* = l/3TF(d 3 2 F)*+4/3W 3 4 F) 

+ 3/5W(d* 2 G) *+ 26/15 W 3 2 H) 

-3/5 W(d 23 F) -2/7 IF(rf 21 D)*- 17/5T7(d 23 F)* 

-12/7W(d 21 G)+4;W(d). (29) 

Formulas for all other terms would require an addi- 
tion to the theory of Backer and Goudsmit. The 
parentages of duplicated terms in r/ 4 and/or d? 
would depend on the particular spectrum under 
consideration, whereas in the formulas for the terms 
given above the coefficients were defined by purely 
mathematical considerations. 

3.2. The l n ~ l V Configuration 

The evaluation of W(y;P r ) when the specification 
7 contains the l n ~ l V configuration, can be carried out 
similarly if tables of coefficients of fractional parent- 
age for nonequivalent electrons are available. 18 In 
general, such tables are not available, and the co- 
efficients must be calculated from tables for equiva- 
lent electrons with the help of transformations 
between different coupling schemes given in section 
2 of III. Again a step-by -step procedure is simplest, 
so that only selections P with n— 1 electrons need be 
considered. The parent term in the l n ~ l configura- 
tion again is assumed, for simplicity, to be in an 
exactly specified state. The single selection con- 
taining the equivalent electrons has an average 
energy ~\{l n - l ySL)V 'S' 'L' ;/ n_1 ) , which is the same as 
the observed energy in the ion W(l n ~ l ySL). The 
only new feature is the selection specified by the 
l n ~ 2 V configuration. The coefficients of fractional 
parentage will have the form 

(l n - 2 (aS l L0r(S''L'') ) lS'L'}l n - l (ySL)l'S'L') 

= (l^iaS^lSLjl*- 1 ySL) 

(S 1 L 1 l(SL)l f S'L'\S 1 L 1 l'(S"L f ')lS'L'). (30) 

On the right, the first factor is the coefficient of 
fractional parentage for equivalent electrons; the 



18 Tables for p 2 s, p 3 s, and d 2 p have been given by S. Meshkov in Phys. Rev. 
91, 871 (1953). With his convention for the normalization of the coefficients, the 
sums over P required in the formulas given here would be replaced by factors 

( ) ; the different selections P are then accounted for during the calculation of 

the tables. To allow what seems to be the most direct use of Racah's tables, the 
coefficients of fractional parentage have been normalized to unity for each selec- 
tion in the present paper. 
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second factor is a transformation coefficient evaluated 
in (5) of III, 

(S 1 L 1 l(SL)l'S'L'\S l L 1 l'(S''L'')lS'L') 

= [(2S+l)(2S" + l)(2L+l)(2L" + l)¥ 



W(S% %S")S'S l )W(Ll'lL")L'L l ). 



(31) 



The W-functions in this expression are calculated 
from (36) of II; some values of the W-f unction have 
also been tabulated. 19 

The particular case where V is an ^-electron is the 
only one we consider in detail. In this case, the 
result of applying (30) and (31) reduces to a fairly 
simple formulation. To evaluate the average 
energies for a term V l ~ l (aSL)sS'L: 

(a) Evaluate W{l n ~ l ySL^ n - 2 ) in terms of the 
observed energies W(l n ~ 2 aSiLi), as already described. 

(b) If the spin Si of the ion term in l"~ 2 does not 
have the same value S' as the original term under 
consideration in Z w-1 $, then add an s-electron to this 
ion term, making the resultant spin S" equal to S, 
the spin of the parent term in l n ~ l ; i. e., change 
W(l n ~ 2 oSiLO to W(l n ~ 2 (aSMsSLJ. 

(c) If the spin Si has the value S", then add the 
s-electron to the ion, and split the latter into two 
terms, each having one of the two spins that are 
allowed, i. e., give $" the two values S and (2$' — S). 
The part having the resultant spin S is given a 
coefficient (2S" + 1)~ 2 , whereas the other part' has a 
coefficient of unitv minus this; i. e., change W(l n ~ 2 
aS,L x ) into [(2S' + 1)" 2 WiP-KaS'LJsS'LO + tt- 
(2S / +l)- 2 ]W(l n - 2 (aS / L l )s(2S'-S)L l )]. 

This formulation could also be used when /' ^0 to 
take care of the spin. Another subdivision would 
then be made among the /^-values, using only that 
part of (31) that depends on the orbital momentum. 
As an example, the second approximation (18) will 
be evaluated for a r/ 3 ( 4 F)s 5 F term. Step (a) is 
already carried out in (25), so that according to 
step (b) 

W((P ( 4 F)s 5 F;d 2 s) 

= l/5W(d 2 (*F)s 4 P)+4/5TF(tf 2 ( 3 *> *F). (32) 

The procedure is repeated for the terms on the right, 
and the results obtained are substituted back in (32) 
to get the average energy for the ds selection: 



W(d\*F)s 5 F;ds) = W(ds 3 D). 



(33) 



This result could also be written by inspection. 
The average energies (32), (33), and (25) are inserted 
into (18), together with the trivial values for 
W(y/P), W(y;d), and W(y;s) ; the sums over P intro- 



19 L. C. Biedenharn, Oak Ridge National Laboratory Report ORNL-109S 
(unpublished). 



duce factors ( ) when P T is the l r configuration 

and ( l ) when P r is the f -1 * configuration: 

W 2 (d\'F)s 5 F) = W(d* *F)+3/5W(d 2 s 4 P) 

+ 12/5 W(d 2 s 4 F) -3/5 W(d 2 3 P) 

-12/5W(d 2 *F)-3W(ds 3 D) + 3W(d) + W(s). (34) 

The following formulas are obtained similarly: 

W 2 (d\±F)s 5 P) = W(d 3 ±F) + 8/5W(d 2 s 4 P) 

+ 7/5 W{d 2 s 4 F)-8/5TF(d 2 3 P)-7 /5W(d 2 3 F) 

-3W(ds *D)+SW(d) + W(s), (35) 

W 2 (d?( 2 R)s l U) = W(d d 2 H)+3/2W(d 2 s 2 F)* 

+3/2 W(d 2 s 2 G)*-Z/2W(d 2 3 F)-3/2TF(rf 2 ! G) 

-3/2W(ds *D)-3l2W(ds ] D)*, 

+ ZW(d) + W(s), (36) 

W 2 (d\ 2 R)s 3 H) *= W(d 3 2 H) +1/6 W{d 2 s 2 F) * 

+4/3W(d 2 s 4 F) + 3/2lF(r/ 2 s 2 G)*-3/2W(d 2 3 F) 

-3/2 W(d 2 l G)-l/2W(ds l T))*-bl2W(ds 3 D) 

+ 3W(d) + W(s). (37) 

Formulas (34) and (35) are valid even when con- 
figuration interaction is present. 

3.3. The l n ~ 2 1' 2 Configuration 

The final configuration considered is l n ~ 2 V 2 . Both 
the l n ~ 2 and the l' 2 configuration are assumed to have 
SL- values that are exactly specified. For the selec- 
tion P n -i that contains only one of the /'-electrons, 
the coefficient of fractional parentage is simply a 
transformation coefficient given by (6) of III : 

(l n - 2 ( y SL)rS 1 L u l / S // UT~\ySL)l /2 (S , L f )S ,f L ,f ) 

= {SLr~{S'L')S''L n \SU'{S l L l )VS''L'') 

= [(2S f + l)(2S 1 + l)(2L / + l)(2L 1 + l)]i 

W(Sh S"\$iS')W(U'L"l'\L l L'). (38) 

This leads to a simple formulation for the l n ~ 2 s 2 con- 
figuration: (d) The parentage of an s-electron in a 
term l n ~ 2 s 2 ySL is divided in the ratio of the multi- 
plicities between the possible parent terms, generally 
two. 
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For the selection P n ~\ that contains both V elec- 
trons, the coefficient of fractional parentage becomes, 

a w - 3 (cvS 1 L 1 )r 2 (S , iO&i2W /, i / l^- 2 (7SL)Z ,2 (S , L , )S ,, I // ) 

^^(aSJ^lSLjl^ySL) 

(S 1 L 1 l(SL)S'L\S''L' / )\S 1 L l S'L'(S 2 L 2 )lS''L''). (39) 

The transformation coefficient is given by (5) of III 
and has a form similar to (31). In a term l n ~ 2 s 2 ySL, 
however, the closed shell of two ^-electrons can be 
disregarded for this selection, so that (39) reduces to 
the procedure already given for a configuration of 
equivalent electrons. 

To evaluate the average energy (6) for a d 2 s 2 3 F 
term the formulation (d) is used: 

W(d 2 s 2 3 F;d 2 s) = 2/SW(d 2 s 4 F) + 1/3 W(d 2 s 2 F), (40) 

W(d 2 s 2 3 F ;ds 2 ) = W(ds 2 2 D) . (41) 

Relation (41) is obtained by inspection; using (d) 
again with (41) leads to the average energy for a 
choice ds, 

W(d 2 s 2 *F;ds)=3/4W(d$ 3 D) + 1/4 W(ds 1 D). (42) 

The result (42) could be obtained also from (40) by 
use of formulations (b) and (c), 

W(d 2 s*F',ds) = W(ds B ~D), 

W(d 2 s 2 F;ds) = l/4W(ds 3 D) +3/4 W(ds X D) . (43) 

When these values are substituted for the observed 
energies W on the right of (40) the result (42) is 
obtained; this checks the consistency of the proce- 
dures and illustrates the use of (c). The values (40), 
(41 ) , and (42) , together with trivial values for W(y ;d 2 ) , 
W(y;s 2 ), W(y)d), and W(y;s) are then inserted into 

(18). The number of selections P r of typeP" is ( ); 

l T ~ l s can be selected in 2 ( ) wavs and l r ~ 2 s 2 

\r—l/ 

(n-2\ 

^_ 2 Jways: 

W 2 (d 2 s 2 3 F) *=4/3TI T (rf 2 s 4 F) +2/3 W(d 2 s 2 F) * 
+ 2W(ds 22 D)*-3W(ds*F>) 
-W(ds 1 B)*-W(d 23 F) 
-W(s 21 S)*+2W(d)+2W(s). (44) 
The following formula is obtained in the same way : 
W 2 (d 2 s 2 3 P)* = 4/3W T (rf 2 8 4 P) +2/3W(d 2 s 2 P)* 
+2 WW 2 D)*-3W(ds 3 D) 

-W(s 21 S)* + 2W(d)+2W(s). (45) 



m 



The formulas for the d 2 s 2 configuration always con- 
tain terms that are likely to be perturbed by con- 
figuration interaction. The s 2 *S term is usually un- 
known, and the applicability of the formulas is 
limited by this also. 

4. Comparison With Experiment 

In all the formulas derived above, the observed 
energies, W, are referred to the core with all four 
electrons removed. Because ionization potentials 
are often not accurate, these absolute energies may 
not be well known. To check the theory of Bacher 
and Goudsmit only differences between terms have 
been calculated. Ionization potentials cancel in the 
differences, and energies W can then be taken as the 
relative-energy levels with respect to the lowest 
term as zero; these are given in a Bureau Circular. 20 
In particular, the following two differences would be 
unaffected by configuration interaction: 

D x = W 2 (d? s 5 F)-W 2 (d* s 5 F) = W(d z *F)-W(d* 4 F) 



+ W(d 2 s 4 P) - W(d 2 s 4 F) - W(d 2 3 P) + W(d 2 3 F) 



(46) 



D 2 = W 2 (d* 5 D)-W 2 (d 3 s 5 F) = 6/5W(d* 4 P) 

+ 9/5W(d* *F)-3/5W(d 2 s 4 P)-12/5^(rf 2 s 4 F) 
-6/5W(d 2 3 P)-9/5 W 2 *F)+3W(ds 3 D) 



+ W(d)-W(s). 



(47) 



In table 1 the calculated and observed values of 
A and D 2 are compared. In the spectra of titanium 
and zirconium sufficient data are available to com- 
pare the Z?i-values. The calculated and observed 
values are in close agreement, the calculated value 
exceeding the observed value by 357 K in titanium 
and by 293 K in zirconium. Separations of terms in 
the same configuration are thus given closely by 
the theory. By assuming that A will be given with 
the same error in more highly ionized spectra, it is 
possible to calculate positions for some terms that 
are not known experimentally, and use the values 
obtained to calculate D 2 in these spectra. 21 Accord- 
ing to table 1 the difference between the two terms in 
different configurations is predicted very badly by 
(47) for titanium, the error in D 2 being 6068 K. 
However, the error is much less for the spectra of 
the higher ionization of vanadium (error 1841 K), 
and it is still less in chromium (error 635 K). The 
error is also less in the homologous sequences of 
zirconium (error 2296 K), where the nuclear charge is 
much greater (Z=40) than in titanium (Z=22), 
but the degree of ionization is the same. The 

20 C. E. Moore, Atomic energy levels, NBS Circular 467, vol. I (1949), vol. II 
(1952). 

21 The d 2 ( 3 P)s 4 P term is unknown in V in and Cr iv; the calculated positions of 
these terms are given in table 1 . If the errors in D\ were to decrease rapidly as 
the errors in D 2 do, then these calculated positions would be about 357 K too high, 
and the calculated values of Di would be too small by about three-fifths of this 
error, or 220 K. This is a small error compared to the total errors found in the 
calculated values of Z>2. However, it is noted later that the errors in D\ probably 
will not decrease. 
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Table 1. Term separations in spectra with'd- and s-electrons 
D x =W(d* s sp)- W(& s 3F) 
D 2 =W(d^B)-W(d 3 8 »F) 



Term 


Ti 


V 


Cr 


Zr 


Observed values 


d<5jr> 

rf3 S 5p 
d3 S 5p 


28881 
6721 
14055 


206 
2926 
13647 


350 
50020 
63261 


22138 
5460 
11103 


£>2 


7334 
22160 


10721 
-2720 


13241 
-49670 


5643 

16678 


d3 4p 

d3 4F 

d2«4p 

d°-s*F 


9452 

1085 

9968 

225 


11668 

336 

(57014) 

44343 


14324 

554 

(119671) 

104637 


9824 
3273 

7860 
758 


rf2 3p 
rf2 3F 

d«3D 


10661 

242 

38277 


13344 

419 

96547 


15868 

662 

167715 


8581 

864 

19063 


d 2 D 

« 2 S 

Di 

Bi 

Di(calc)-7)i(obs) 
7)2(calc)-D 2 (obs) 


230 
80379 


372 

148100 


574 
227775 


750 
38258 


Calculated differences 


7691 

28228 


(11078) 

-879 


(13598) 
-49035 


5936 
18974 


Errors 


357 
6068 


« (357) 
1841 


a (357) 

635 


293 
2296 



■Values in parentheses are calculated to make the error in D\ the same for all 
pectra with 3d- and 4s-elcctrons. 



improvement from titanium to vanadium to chro- 
mium is especially remarkable in view of the fact that 
the errors in D 2 form a rapidly decreasing geometric 
sequence, whereas D 2 itself starts at a large positive 
value and decreases to a twice as large negative value. 

Racah 22 has pointed out that the energies would 
be expected to vary more smoothly if the degree of 
ionization (rather than the atomic number) were the 
same in the different spectra. The results obtained 
by use of this empirical modification of the Bacher 
and Goudsmit theory are given in table 2. Except 
for the first spectra, the errors in calculated values 
of D 2 are on the order of 100 K; part of the good 
agreement must be fortuitous, however. 

A partial explanation of the poor agreement ob- 
tained in the first spectra may be that the d? 4 P (of 
Sc r and Y i) is perturbed by interaction with the 
term dsd'f 4 F (in the spectra with 4rf- and 5s-electrons, 
this has the added effect of making the calculated 
position for the d 2 3 F too high, as seen from the small 
separation of the d 2 3 P and d 2 3 F in Sr i). This effect 
can be estimated in the spectra containing 3d- and 
4s-electrons by calculating the position of d? 4 P so 
that the error in T)\ will be 170 K as it is in the second 
spectra. When this calculated position (40048 K) 
is used, the error in D 2 is still very large (—5769 K). 
Racah has rioted that the observed position of d 2 D 
in K i is about 5117 K lower than expected from the 
behavior of the difference between baricenters of the 

22 G. Racah, Bui. Research Council Israel III, 290 (1954). I am indented to 
Professor Racah for making a preliminary copy of this manuscript availahle and 
for discussing the relationship of his results to those of this paper. See also M. 
A. Catalan, F. Rohrlich, and A. G. Shenstone, Proc. Roy. Soc. [A] 231,421 (1954). 



Table 2. Term separations in spectra with d- and s-electrons 



Term 


Spectrum 


I 


II 


III 


IV 


I 


II 


Observed values 


d*5D 

rfs * *F 
rf3 S 5P 

D x 
Do, 

#<p 

rf3 4F 
rf2 S 4p 
rf2 S 4F 

rf 2 3P 
rf2 3F 
rf«3D 

S2S 

D x 
D 2 

Di(calc)-ZMobs) 
ZMcalc)-Z> 2 (obs) 


(Ti) 28881 
6721 
14055 

7334 

222160 

(Sc) 36540 
33846 
17969 
11610 

(Ca) 48551 
43494 
20357 

(K) 21535 



(V) 206 
2926 
13647 

10721 
-2720 

(Ti) 9452 

1085 

9968 

225 

(Sc) 12128 

4909 

105 

(Ca) 13687 



(Cr) 350 
50020 
63251 

13241 

-49670 

(V) 11668 

336 

(56841) 

44343 

(Ti) 10661 

242 

38277 

(Sc) 118 
25537 


(Mn) 537 
L 12296 


(Zr) 22138 
5460 
11103 

5643 

16678 

(Y) 32210 
29605 
15385 
11277 

(Sr) 44662 

(43989) 
18254 

(Rb) 19355 



(Nb) 772 
3307 
11034 

7637 
-2625 

(Zr) 9824 
3273 

7860 

758 

(Y) 14048 
8429 
1193 

(Sr) 14724 





-111759 

(Cr) 14324 

554 

a 119671 

104637 

(V) 13344 
419 

96547 

(Ti) 230 

80379 


Calculated differences 


3996 
12181 


10891 
-2613 


(13411) 
-49738 




(6040) 
-3012 


8034 
-2582 


-112020 


Errors 


-3338 
-9979 


170 
107 


b (170) 
-68 




b (397) 
-19690 


397 
43 


-261 



a From table 1 calculation. 

b Values in parentheses calculated to make error in A have the indicated value. 



301746—54- 
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d n ~ l s and d n configurations in Ca i, Sc i, Ti i, and 
V i, and this may explain the poor agreement ob- 
tained in first spectra. 

5. Relations Between Radial Parameters 

All formulas of the Bacher and Goudsmit theory 
can be used also to calculate the matrix elements of 
scalar operators; this has been noted in specific cases 
in the discussions of formulas (21) and (22), and it 
was also pointed out in BG. If the term values could 
be expressed in terms of such matrix elements, and 
if the zero-order functions were the same for all 
configurations, then the formulas derived from (16) 
to (19) (such as (46) and (47)) would be identities if 
the W's were replaced by matrix elements of these 
operators. 

The matrix elements of the operators are linear 
combinations of a limited number of radial integrals 
For the Slater theory, these combinations are well 
known (most of the formulas necessary for the pres- 
ent paper are contained in II). The observed sepa- 
rations of terms in the same configuration are pre- 
dicted well by the Slater theory (i. e., with an accu- 
racy equal to that obtainable with the Bacher and 
Goudsmit theory) if an empirical correction pro- 
portional to Z(L+1) 23 is added to account for the 
polarization energy (see footnotes 8 and 9). .Radial 
integrals defined by zero-order one-electron func- 
tions having the same n- and /-values must be con- 
sidered independent of each other when in different 
configurations to get this good agreement, however; 
if the Bacher and Goudsmit theory were directly 
applicable, all such integrals would have to have the 
same value. The derivation of the theory, given in 
section 2, shows how much of the observed variation 
of these radial integrals will be accounted for by 
the Bacher and Goudsmit theory in any given 
approximation. 

To illustrate, formulas for the W's that enter into 
(46) are given here according to the notation used 
in section 2: 

W(d 3 s 5 P) - A(d 3 s) —3G 2 (fP s) -I6a(d 3 s) 

W(d* s 5 ¥)=A(d 3 s)-!5B(d* s)-3G 2 (d 3 s)-6a(d* s) 

W(d 3 *F)=A(d 3 )-Wa(d 3 ) 

W(d 3 *F)=A(d 3 )-l5B(d 3 )-6a(d 3 ) 

W(d 2 s *F)=A(d 2 s)+7B(d 2 s)-2G 2 (d 2 s)-10a(d 2 s) 

W(d 2 s*Y) = A (d 2 s) - SB (d 2 s) -2G 2 (d 2 s) 

W(d 2 *F)=A(d 2 )+7B(d 2 )-10a(d 2 ) 

W(d 2 3 F)= A(d 2 )-SB(d 2 ) (48) 

23 The L(L-\-l) correction is represented as the matrix element cf the scalar 
n 
operator 2« 2 h'U- The value of this for a term in the l n s k configuration 

is a[L(L+l) — nl(l-\-l)]. In the present work, this, rather, than simply «L(L+1), 
is added to each term; this change in the form of the correction is compensated by 
a change in the value of the constant parameter A. The modification has been 
made so that the relationships obtained later between the ^4-values in different 
configurations will have a simpler form. 



If these values are substituted in (46), then the latter 
is an identity if the B's and as in different configura- 
tions are assumed equal. Relation (46) will still be 
satisfied if the following two less restrictive conditions 
hold: 

B{d? s )=B(cP) +B(d 2 s) -B(d 2 ) , 



a(d z s) = a(d*)+a(d 2 s)-a(d 2 ). 



(49) 



It is sufficient for B—2/3a to satisfy a relation of 
similar form for (46) to be satisfied, but if separations 
of all pairs of terms in the d 3 s configuration are to be 
correct, then both relations must hold and in addition 
the following two : 



C(d* s) = C(d?) + C(d 2 s) - C(d 2 ) , 
G 2 (d 3 s)=2G 2 (d 2 s)-G 2 (ds). 



(50) 



Formulas of this type can be derived easily by 
inspection, and there is no need to consider fractional 
parentage. To obtain the second equation in (50), 
for instance, note that in (18) the combination of a 
d- and an ^-electron (needed to define G 2 ) occurs three 
ways in W 2 (d z s) ; it occurs two ways in each of three 
possible selections P 3 that contain the s-electron (in 
the d 2 s configuration) ; it occurs only once in each of 
the three possible selections P 2 that contain the 
s-eiectron (in the ds configuration). As far as the 
relation for the parameter G 2 is concerned, this direct 
interpretation of (18) leads immediately to (50) 
(multiplied by a factor of three). As B, C, and a 
are all defined by the combination of two rf-electrons, 
a similar derivation leads to the same form of equa- 
tion for all three of these parameters, and this 
equation is the one represented by (49) and the first 
equation of (50). In this way it can be shown that 
the Bacher and Goudsmit second approximation 
(relation (18) with 72 = 4) will predict correctly the 
separations of terms in the same configuration if the 
following relations hold between the parameters : 

B(d i ) = 2B(d 3 )-B(d 2 ) (51) 

B(d z s)=B{d?)+B{d 2 s)-B(d 2 ) (52) 

B(d 2 s 2 )=2B(d 2 s)-B(d 2 ). (53) 

Three similar relations must be satisfied by the C's 
and three more by the a's. Also, it is necessary 
that 



G 2 (d 3 s)=2G 2 (d 2 s)-G 2 (ds). 



(54) 



The derivation from (18) fails to bring out the fact 
that the two-electron parameters must have a linear 
variation if the formulas are to be valid. In general, a 
two-electron parameter defined between electrons 
with orbital angular momentum I and V in the 
configuration y will be specified by (y)G(ll f ). Ac- 
cording to the derivation with the second series of (9) , 
only a first-power variation with the number of 
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electrons in 7 can be accounted for by the Backer 
and Goudsmit second approximation. In particular 
then, any one of the parameters B, C, or a must 
satisfy a relation of the form 



(d r ~ k s k )G(dd)=A+rB+W. 



(55) 



It is easy to verify that (51) to (53) are satisfied 
when a formula of this kind is applicable. 

It is possible to generalize (55) to include all 
similar spectra rather than just those of the ions 
of the same atom, in two-, three- and four-electron 
spectra. If q is the degree of ionization, this general- 
ization will have the form (d r ~ k s k )G(dd)=A J r rB+ 
kC+qD. The difference A defined by (46) (and 
equal to 15/?— 10a according to the linear theory) 
should also be expressible by a formula of this type. 
Racah 22 has fitted 34 experiment ally observed P-F 
differences in spectra with 3d- and 4s-electrons with 
a mean error of 340 K by use of the formula 

A = 2679+853r4-1334t + 3158g-lllr/. (56) 

This formula is adjusted to spectra with eight 
different ionizations, and without a quadratic term 
the mean error would be increased to about 520 K. 
According to this formula, the differences D\ cal- 
culated from (46) will be too great on the average 
by about 220 K. Modifying the theory so that 
spectra with the same ionization are utilized tends 
to eliminate this systematic error. This could also 
be accomplished by using the \\\ approximation. 

The original series (9) will often have to be satisfied, 
for a given y, by more than one value of a, and this 
leads to the conclusion that the ratios of similarly 
defined radial parameters in different spectra of a 
given element should be equal. In particular, the 
following ratio must be constant in configurations 
with d- and s-electrons: 



B(d r s k ) 



C(d r s k ) 



^constant. 



(57) 



This is not really a necessary consequence of the 
theory developed in section 2, however, as the series 
(9) can be reinterpreted in terms of radial parameters 
that independently satisfy relations such as (55). 
When irrational roots are considered, it is shown in 
section 7 that (57) is a necessary condition for con- 
sistency of this interpretation, in terms of a linear 
theory, with a direct application of the method of BG. 
In the same way that the differences such as (46) 
were replaced by the formulas (51) to (54), the 
differences similar to (47) can be replaced by rela- 
tionships between the A values. The Backer and 
Goudsmit second approximation will give the sepa- 
rations of terms in different configurations correctly 
if the following two relationships between the ^4's 
are satisfied: 

A(d') -A(d" s) = 3A(d*)-3A(d 2 s) -3A(d 2 ) 

+ 3A(ds)+A(d)-A(s) (58) 



A{d z s) -A(d 2 s 2 ) =A(d 3 ) +A(d 2 s)-2A(ds 2 ) 
-2A(d 2 )+A(ds) +A(s 2 ) +A(d) 

-A(s)+2G 2 (d 2 s)-2G 2 (ds) (59) 

By redefining the ^1-values suitably, 24 the dependence 
on G 2 may be omitted from (59); equations (58) and 
(59) may therefore be regarded as two relations 



between the differences A (d r ~ k s k ) — A (d r ~ 



l ). In 



the second approximation, the parameter A must 
satisfy a cubic relation, and these differences then 
satisfy quadratic equations having the general form 



A(d r ~ k s k )-A(d 



r-k-l Q k+l\ 



)=D+Ek + Fr+Gkr+Hr 2 . 

(60) 



The k 2 dependence is omitted because it adds nothing 
extra when k has only the two values and 1. It is 
easy to verify that (58) and (59) are satisfied when 
a formula like this is applicable. This formula can 
be generalized to apply to all similar spectra in a 
manner similar to that used for (55) if the "A" is 
suitably defined in the spectra with more than four 
electrons. If third (or fourth) power terms were 
important, a constant (or linear) error would result. 
The errors of I) 2 in table 1 change too rapidly to be 
explained in this way, however. 

In the one-electron, two-electron, and three- 
electron spectra considered in section 4, there are 
five experimentally determinable differences between 
yl-values, and the five unknown constants in (60) 
can be evaluated from these observed values. Simi- 
larly, the three constants in (55) can be evaluated 
from observed separations of terms in the two- 
electron and three-electron spectra. The Bacher 
and Goudsmit second approximation is then equiva- 
lent to using the formulas obtained in this way to 
extrapolate values in the four-electron system. 

This equivalence has been verified by using all 
available experimental data in titanium and zir- 
conium to evaluate the radial parameters for all 
spectra with least squares. 25 The results confirmed 



24 The parameter G 2 appears in (59) because an -s 2 subshell is neglected in defining 

the - -1 -values; thus the formula for a d 2 s 2 Z F term would have the same form as 
the last relation in (48). If the effects of the ^-electrons in the * 2 subshell were 
explicitly expressed, as done when only a single 8-electron Is present, then 

2G 2 (d 2 s 2 ) would be subtracted from the formulas for all terms in the d- s 2 con- 
figuration. This has not been done since the effect would be included in . i(d 2 .s 2 ) 
if the parameters were evaluated by least, squares and G?(d 2 s 2 ) could not be 
evaluated. This applies also to G 2 (ds 2 ). The G's that cannot be evaluated 
satisfy the relation c7 2 (f/ 2 s 2 )— G 2 (ds 2 ) = G 2 (d 2 s)-G 2 (ds) = G 2 (d 3 s)-G 2 (d 2 s). 

25 The a *F term of Tit and the a *S and c »F terms of Zr I were excluded. The 
a U? term is probably not real. It was predicted 2000 K lower than observed when 
it was not included in the least-squares calculation. As long as it was included, it 
and several other terms were in error by 500 to 600 K, regardless of uncertainties 
in effects of configuration interaction. It is not well-established experimentally, 
as only three weak lines were observed; one strong line used to identify it is a 
blend (H. N. Russell, Astrophys. .1. 66, 283 (1927)). The a is term was calculated 
8000 K higher than observed; as *S terms arc always hard to establish, it is fairly 
certain that this term is also unreal. The c 5 F term was calculated 5000 K higher 
than observed; although highly unlikely, it is not certain that consideration of 
of configuration interaction could not explain the result. However only one strong 
transition has been observed from each level of this triplet; five other confirming 
transitions are all weak (C. C. Kiess and II. K. Kiess, BS J. Research 6, 621 (1931) 
RP296). The a 5 D, b *Q, and d sp terms in the rf 4 configuration of Ti I are observed 
much closer together than called for by the theory; this may be because cf con- 
figuration interaction, or it may represent a failure of the W 2 approximation. 
In addition to the experimental values listed in Atomic Energy Levels, the re- 
cently located terms d 2 D of Zr n and the a ! S and b 'S terms of Zr in (C. O. Kiess, 
J. Opt. Soc. Am. 43, 1024 (1953) )were also available. I am indebted to Dr. Kies 
for furnishing these data before publication. 
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the conclusion already obtained from the limited 
amount of data utilized in table 1. In all spectra 
of both elements, the terms could be fitted with 
absolute errors less than 400 K. The parameters 
obtained could be adjusted to fit the relations con- 
tained in (51) to (57) without significantly changing 
the original agreement. This verifies the conclusion 
that the separations of terms in the same configura- 
tion are given by the Bacher and Goudsmit method 
with an absolute error less than 400 K in titanium and 
zirconium. The discrepancy in calculating A(d*) — 
A(<P s) from relation (58) was found to be 6560 K 
in titanium and 2086 K in zirconium. These values 
are approximately the same as the errors found for 
D 2 in table 1, which verifies the conclusion that 
calculated separations of terms in different configura- 
tions may be considerably in error. 

If the ionization potentials are accuratley known 
the aboslute values of the A's can be calculated 
(rather than just the differences (58) and (59)), 



A(d*)=4A((P)-6A(d 2 )+4A(d)+ 8 



(61) 



A(d 3 s) =A(d 3 )+ 3A(d 2 s) -3A(d 2 ) 

-ZA(ds)+ A($)+3A(d)+ 8 (62) 
A(d 2 s 2 ) = 2A(d 2 s)+ 2A(ds 2 ) -A(d 2 ) -4A(ds) -A(s 2 ) 



+ 2A(d)+ 2A(s)-2G 2 (d 2 s)+2G 2 (ds) + 8, (63) 



where 



5 = 71-3/*+ 3/3-/4. 



The ionization potentials appear only in the combi- 
nation 8; /1 is the principal ionization potention of 
the four-electron spectrum (T\ 1, for instance); I 2 
is that for the three-electron spectrum (as T'\ 11), etc. 
If (61) to (63) are utilized, errors in the theory will 
not cancel as they may in the differences (58) and 
(59). In titanium the error of (61) in calculating 
A(d^) is found to be 24130 K, and the error of (62) 
in A(d 3 s) is 1757 K; on the other hand, the corre- 
sponding errors in zirconium are — 2100 K and —4186 
K. The inconsistent behavior of the errors in these 
two calculations probably arises from the well-known 
need for better ionization potentials, particularly 
in the one- and two-electron spectra. 

6. Weak Configuration Interaction 

The proof of the method of Bacher and Goudsmit 
is established with second-order perturbation theory, 
so that the theory is not directly applicable if con- 
figuration interaction is too strong to be evaluated 
with perturbation theory. However, it may not be 
directly applicable even if the spectrum has weak 
configuration interaction that can be evaluated with 
perturbation theory. 

For example, consider interaction between terms 
of d n and those of d n ~ l s and assume that second-order 
perturbation theory can be used. The perturbation 
of, say, a d 3 2 G term by a d 2 s 2 G term is evaluated 
in terms of a characteristic interaction parameter 



associated with the d z and d 2 s configurations; this is 
designated H 2 (d 3 —d 2 s) in agreement with the nota- 
tion in (80) of III. For arbitrary n, there is a similar 
parameter defined between d n and d n ~ l s. Partly 
with the help of table XX of III the following two 
perturbations are evaluated: 



A(d 32 G) = 



A(d 2l B)- 



50[H 2 (d 3 -d 2 s)] 2 
W(d 32 G)-W(d 2 s 2 G) 

14:0[H 2 (d 2 -ds)] 2 
: W(d 21 B)-W(ds l D) 



(64) 



(65) 



By considering terms marked with asterisks in (28) , 
it follows that configuration interaction between d n 
and d n ~ l s will be correctly accounted for by the 
Bacher and Goudsmit theory if the following relation 
exists between the two perturbations (64) and (65) : 



6/5A(rf 3 2 G) -3/7 A(d 2 X D) =0. 



(66) 



According to the perturbation calculation used by 
Bacher and Goudsmit, the energy difference W(d 3 
2 G)-W(d 2 s 2 G) would equal the difference W(d 2 
X D) — W(ds l B) and the parameter H 2 (d 3 —d 2 s) would 
equal H 2 (d 2 —dx). If these two conditions held, then 
(66) would clearly be satisfied. Relations similar to 
(66) can be found from other formulas (such as (29) , 
(36), and (37)). By assuming that all energv differ- 
ences W(d n aSL) — W{d n - 1 sa'SL) are equal* and all 
parameters H2(d n —d n ~ l s) are identical, matrix ele- 
ments of configuration interaction can be calculated 
with the Bacher and Goudsmit theory by solving a 
set of simultaneous equations, although such a pro- 
cedure would not be the most practical. Similar 
considerations would apply 7 to any pair of configura- 
tions; in particular, the interaction between d n s and 
d n ~ l s 2 is interesting because the parameter defined 
between these two configurations would be equal to 
the parameter defined between d n and d n ~ l s accord- 
ing to this direct application of the Bacher and Goud- 
smit theory. 

Alternatively, (66) could be satisfied if the para- 
meters H 2 (d 3 —d 2 s) and H 2 (d 2 —ds) accidentally had 
the same ratio as the square roots of the energy 
denominators in (64) and (65). To get consistency 
for all formulas similar to (66) that could be derived, 
it is necessary that a representative average energy 
difference exist that approximates separations of 
terms with the same S- and Z-values in the two 
configurations (in many cases no average value will 
adequately represent the data). An easily defined 
value for this average is the difference in the ^1-values 
of the two configurations. In table 3 calculated 
iZ 2 -values and differences between ^4-values have 
been collected. In a given atom, the differences 
between A- values in this table are greater the 
greater the ionization; H 2 is also larger in the spectra 
of higher ionization, according to the results obtained 
for vanadium. 

An extension of this argument indicates that the 
polarization energy arising from interactions with 
configurations in which only one electron is excited 
may be partially included in the second approxima- 
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Table 3. Average term separations and interaction parameters 



A(s)-A(d) 


Titanium 


Vanadium 


Chromium 


ZircoTiium 


A(a)-A(fi) 


Ih 


A(a)-A<fi) 


Ih 


A(a)-A(fi) 


Ih 


A(a)-A(0) 


Ih 


80150 




147730 




227200 




37508 




A(ds)-A(<V) 
A(s*)-A(ds) 


34470 




92460 




162720 




16700 
20304 




A(dU)-Am 
Hi 
A(ds*)-A(d? 8) 


a -896 
a 16980 


a 172 


a 42102 




103100 




b -1260 
b360O 


t-387 


A(dzs)-A(d*) 

H 2 

A(<P 8*)-A(d* s) 


-17140 
a -13103 




a 6332 


a 153 


52424 




-13230 
-11700 




A(d*s)-A(d*) 
Ih 

A{d*8*)-A{d°>8) 
Ih 






« -21335 

o -8759 


c 106 


o 11875 
30400 


" 150 







a A. Many, Phys. Rev. 70, 511 (1946). 
b C. W. Uflford, Phys. Rev. 44, 732 (1933). 
o A. Schweizer, Phys. Rev. 80, 1080 (1950). 



tion of BG. However, it is unlikely that the inter- 
action parameters would behave exactly as required 
so that it might not be well accounted for. 

It is clear from table 3 that in the titanium spectra 
of moderate ionization, the Bacher and Goudsmit 
method would not account for the configuration 
interaction. Even though the interaction is weak 
enough to be evaluated with second-order perturba- 
tion theory, the calculations discussed in section 5 
show that errors of 300 to 700 K would be introduced 
into formulas (28), (29), (36), and (37). Whenever 
the average term separation is small or negative, 
which is the case in Ti i and Ti n, no set of values for 
the interaction parameters can produce agreement. 
When the configuration interaction is so strong that 
perturbation theory cannot be used, the errors are 
even larger. For zirconium, errors in excess of 1000 
K result when the formulas are applied directly. 
The differences in A- values are all large and positive 
in chromium, so that the configuration interaction 
would probably be well accounted for in the second 
approximation; in this case, however, the interaction 
is weak and could be neglected in the linear theory. 

7. Irrational Roots 

This section deals with one of the simplest exam- 
ples in which irrationality of the energy levels plays a 
part. As a preliminary, the two 2 D terms of the d 3 
configuration are discussed. 

The determination of the energy levels of the two 
2 D terms of d 3 was first carried out by Ufford and 
Shortley. 26 They chose two mutually orthogonal 2 D 
eigenfunctions, which they characterized as a 2 D and 
b 2 D; in the more recent notation of III these same 
two eigenfunctions are designated as 2 iD and 2 3 D, 
respectively. 27 In terms of this choice of eigenfunc- 
tions, the electrostatic energy matrix was evaluated 

26 C. W. Ufford and G. H. Shortley, Phys. Rev. 42, 167 (1932). See alsc, E. U. 
Condon and G. II. Shortley Theory of atomic spectra, p. 233 (Cambridge, 1951). 

27 In III the classification of terms in d* is carried out by introducing the senior- 
ity number; tin- is given as the prefixed subscript in the notation above. 



and found (o Lave the form 

M JD ( A + 7B + 7C 3 V 2T£ \ 

The energy eigenvalues of the two 2 D states of d 3 are 
the characteristic values of this matrix; these values 
(as well as the observed energies for the 2 D terms) 
are specified by IF( 2 D+) and W( 2 D"). 

The matrix S that diagonalizes (67) will have the 
form 

s =(-c 2 c) {Q+cl=1) ' (68) 

The elements of S can be specified by the condition 
for the vanishing of the nondiagonal elements of 
S M S~\ 



[U3(d 3 )+4C(d 3 )]C 1 C 2 -(C 1 2 -C 2 2 )S^21B(d 3 )=0. (70) 



Although ^(iD) and W(aD) have no physical reality, 
they will be used subsequently to indicate the 
energies that would be associated with the diagonal 
elements of the matrix (67) . The following formulas 
are obtained from (69) : 



WQB) = C! W( 2 D-) + Ci W( 2 D+), (71) 

W(ID) = C! W( 2 B + ) + Ci W( 2 D-) . (72) 

If the polarization energy can be accounted for with 
the L(L-\-\) correction, then all the preceding dis- 
cussion will still apply even when polarization energy 
is considered if the constant A in (67) is suitably 
redefined. 

By direct application of the methods in section :\, 
the energy for, say, a d 4 3 G term, can be expressed in 
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the form 

W 2 (d* 3 G)* = 1 0/2 1 W(d 3 2 3 D)* + 3/4 W(d 3 2 F)* 

+ 4/3W(d 3 "F) + 99/140W(d 3 2 G)* 

+ 1 1/1 5 W(d* 2 H)- 11/10 W(d 2 3 P) 

-9/14W(d 21 D)*-29/10W(d 2z F) 

-19/14W 7 (rf 21 G) + 4W^). (73) 

If the linear theory is valid, Ti'XID) in (73) can be 
replaced by the value (71); the approximation 
obtained in this way will be specified by a prime. In 
the exact application of the Bacher and Goudsmit 
theory, (71) can still be used to replace W(d? IT)), 
but C\ and C 2 are not necessarily specified by (70). 
Also, instead of evaluating W(d 3 %D;d 2 ) which was 
done in deriving (73), it is necessary to evaluate 
W(d 3 \C, 2 D~ + C 2 2 B + ];d 2 ). The following is the 
correct formula which is obtained in this way: 

W 2 (d* 3 G)= W 2 (d± 3 G)-e(d 4 3 G) (74) 

e(d* 3 G) = C\ CJ[8/21 W(d 2 ! S)-2/7 W(d 2 3 P) 
-50/147W 7 (rf 21 D) + 2/7W^/ 23 F) 

- 2/4 9 W(d 2 1 G)] - C ^ C ^-°^ [2lW{d 2 3 P) 



-15W(d 2 l D)-2lW(d 2 3 F)+15W(d 2 'G)}. 

(75) 

Replacing the d 2 energies by their Slater formulas 
(including the L(L + 1) correction) leads to the 
following approximate form for (75) : 

e(d* 3 G)=10l2lC 1 C 2 {[4B(d 2 ) + 4C(d 2 )]C 1 C 2 -(C 2 1 

-<7I)3V2lB(d 2 )} (76) 

It follows that when C\ and C 2 are determined from 
(70) then (74) will be equivalent to (73) if the ratio 
of B to C is the same in d 3 as it is in d 2 . To establish 
the general correctness of the method of applying 
the Bacher and Goudsmit theory which was devel- 
oped in section 5, it is necessary that in general the 
ratios of similarly defined radial parameters be the 
same in all spectra of a given atom. 

It is well known that the ratio of B to C is very 
nearly constant. Values of this ratio, obtained 
mainly from the calculations discussed in section 5, 
are given in table 4. Because of the large errors that 
may be present in the theory (±400 K), the values 
of r 2 , obtained by fitting the observed data exactly, 
are not well defined; to a lesser extent the d 3 and d 3 d 
configurations also have too few terms to define n 
precisely. The data as a whole (giving strong weight 
to the two well determined r 5 -values) support the 
conclusion that r n is a slowly decreasing function of 
?i. The exception apparent in titanium is less 
noticeable if the irrational roots are omitted from 
the least-squares calculations for 3d 3 and 3d 3 4s 
parameters (the error in fitting the irrational terms 
is then excessive, however). 

To estimate roughly the effect that could result 
from inequality of this ratio, (76) can be put in the 



following form with the help of (70), and by use of an 
approximate value of r equal to 0.29: 



€ (d 43 G) = 40/21 



e?) 



C\ Ct C(d 2 ) 



= l.lQ(r 3 -r 2 )C(d 2 ). 



(77) 



In Ti in, C(d 2 )=2505 so that if r 3 — r 2 =0.05 then 
e=145 K. This is likely to overestimate the max- 
imum error that can result, yet within the accuracy 
of the theory it probably could be neglected. How- 
ever, the same treatment for the r/ 4 X F term shows 
that e will have exactly 3 times the value (75), and 
for the rf 4 3 D term the error will be 2.4 times this 
value, so it is not certain that this effect will always 
be negligible. The ratio 6 2 /C need not be constant 
(it varies considerably, from 0.4 to 0.8 usually). 

Table 4. BjC ratios 

r n is the ratio of B(d n ) to C(d n ), R(d» s) to C(d n s), or P(d» 8*) to C(d n s 2 ), accord 
ing to whichever definition has meaning for the spectrum specified. 



Ti 1 


n 


rz 


Th 


0.289 


«0. 343 (0.319) 




Ti ii 


.290 


.317 (0. 290) 




Ti in 


.291 






Fe in a 






0.271 


Zr i 


. 304 


.279 




Zr ii 


.289 


.287 




Zr in 


.297 






Mo i b 






. 257 



a See footnote 8, page 38. b See footnote 3, paee 35. 

Values in parentheses are calculated with irrational 
roots omitted. 

Though the method of BG introduces new non- 
linear effects if the parameter ratios are not equal, 
these effects would probably represent errors in the 
method. Because of this, and because the ratios 
are so nearly equal, it is not expected that a correc- 
tion with this variation will produce improved 
agreement when added to the linear theory. How- 
ever, irrational terms often cannot be fitted by the 
linear theory with the accuracy that is obtainable 
for rational terms 28 and some additional correction 
of this form may be needed. This correction would 
not change the agreement obtainable for the rational 
terms of c/ 4 , but it would introduce corrections into 
the irrational terms. These considerations may be 
less academic for spectra with p- and ^-electrons 
where there are more parameter ratios that may 
differ, and where most of the terms are irrational. 29 



I am indebted to C. C. Kiess for suggesting the 
problem on which this work is based, and I thank 
also G. Racah and C. W. Ufford for their comments 
on the original manuscript of this paper. 

28 R. E. Trees, Phys. Rev. 85, 382 (1952). 

29 Meshkov has introduced a method which combines the Slater theory with 
a Bacher and Goudsmit first approximation, and has obtained improved agree- 
ment for the 3d 2 4p configuration of Ti u (footnote 18). A simple extension of 
the linear theory to include correlations between p- and d-electrons may explain 
this improvement, but the possibility that new types of crorections are needed 
to explain it, cannot yet be excluded." His theory also leads to better agreement 
than is obtainable with the Slater theory alone, in the d 3 configuration of V in 
and the r/ 4 configuration of V n (Phys. Rev. 93, 270 (1954)). In the d* configura 
tion his final mean deviation of ±663 K is still rather large; this may be because 
there are a large number of irrational terms present in this configuration or 
because configuration interaction was neglected. See also S. Meshkov and C. 
W. Ufford, Phys. Rev. 94, 75 (1954). 

Washington, January 8, 1954. 
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